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TECBNICAL  NOTE  21k0 

THEORETICAL  ANTISYMMETRIC  SPAN  LOADING  FOR  WINGS  OF 
ARBITRARY  PLAN  FORM  AT  SUBSONIC  SPEEDS 
By  John  DeYoung 


SUMMARY 


A  simplified  lifting-surface  theory  that  includes  effects  of  com- 
pressihility  and  spanwise  variation  of  section  lift— curve  slope  is  used 
to  provide  charts  with  which  antisymmetric  loading  due  to  arbitrary  anti¬ 
symmetric  6Lngle  of  attack  can  be  found  for  wings  having  symmetric  plan 
forms  with  a  constant  spanwise  sweep  angle  of  the  quarter-chord  line. 
Consideration  is  given  to  the  flexible  wing  in  roll.  Aerodynamic  char¬ 
acteristics  due  to  rolling,  deflected  ailerons,  and  sideslip  of  wings 
with  dihedral  are  considered.  Solutions  are  presented  for  straight- 
tapered  wings  for  a  range  of  swept  plan  forms. 


INTRODUCTION 


Reference  1  has  been  for  many  years  the  standard  reference  for  esti¬ 
mating  the  stability  and  control  characteristics  of  wings.  The  lifting¬ 
line  theory  on  which  this  work  was  based  gave  generally  satisfactory 
results  for  straight  wings  having  the  aspect  ratios  considered;  however, 
the  use  of  wing  sweep  combined  with  low  aspect  ratio  has  made  an  exten¬ 
sion  of  this  work  desirable.  Lifting— line  theory  cannot  adequately 
account  for  the  increased  induction  effects  due  to  sweep  and  low  aspect 
ratio;  consequently,  it  has  been  found  necessary  to  turn  to  the  more 
complex  lifting— s-ur face  theories. 

Of  the  many  possible  procedures,  a  simplified  lifting-surface  theory 
proposed  by  Weissinger  and  further  developed  and  extended  in  reference  2 
has  been  found  especially  suited  to  the  rapid  computation  of  character¬ 
istics  of  wings  of  arbitrary  plan  form.  Comparisons  with  experiment  have 
generally  verified  the  theoretical  predictions.  In  references  3  and  4, 
this  method  has  been  used  to  compute  for  plain,  unflapped  wings, the  aero¬ 
dynamic  characteristics  dependent  on  symmetric  loading.  The  same  simpli¬ 
fied  lifting— surface  theory  can  be  extended  to  predict  the  span  loading 


2 


NACA  TN  2lk0 


resulting  from  antisymmetric^  distribution  of  the  wing  angle  of  attack. 
From  such  loadings,  the  damping  moment  due  to  rolling,  the  rolling 
moment  due  to  deflected  ailerons,  and  the  rolling  moment  due  to  dihedral 
angle  with  the  wing  in  sideslip  can  be  determined.  A  recent  publication 
(reference  5)  makes  use  of  the  simplified  lifting— surface  theory  to  find 
span-loading  characteristics  of  straight-tapered  swept,  wings  in  roll  and 
loading  due  to  dihedral  angle  with  the  wing  in  sideslip.  Experimental 
checks  of  the  theory  for  the  daii5)ing— in— roll  coefficient  and  rolling 
moment  due  to  sideslip  were  very  favorable.  The  range  of  plan  forms 
considered  in  reference  5  is  somewhat  limited  and  aileron  effectiveness 
was  not  included.  The  loading  due  to  aileron  deflection  normally  involves 
excessive  labor  when  computed  by  means  of  the  si]25)lified  lifting— surface 
theory;  however,  development  of  the  theory,  presented  in  reference  6, 
that  deals  with  flap  and  aileron  effectiveness  for  low-aspect-ratio  wings 
provides  a  means  by  which  the  simplified  lifting-surface  method  can  be 
used  to  obtain  spanwise  loading  due  to  aileron  deflection. 

It  is  the  purpose  of  the  present  analysis  to  provide  simple  methods 
of  finding  antisymmetric  loading  and  the  associated  aerodynamic  coeffi¬ 
cients  and  derivatives  for  wings  with  symmetric  plan  forms  limited  only 
by  a  straight  quarter— chord  line  over  the  semi  span.  Means  will  be  pre¬ 
sented  for  finding  quickly  the  aerodynamic  coefficients  of  span  loading 
due  to  rolling,  of  span  loading  due  to  deflected  ailerons,  and  of  span 
loading  due  to  sideslip  of  wings  with  dihedral.  Flexible  wings,  when 
the  flexure  depends  principally  on  span  loading  as  in  loading  due  to 
rolling,  can  be  included  in  the  analysis. 


NOTATION 


A 

b 


aspect  ratio 


wing  span  measured  perpendicular  to  the  plane  of  symmetry,  feet 


c  ^wing  chord,  feet 

Ca  ^aileron  chord,  feet 

Cav  ^mean  wing  chord  ^  ,  feet 

Cj  local  lift  coefficient  ^ 

iThe  word  antisymmetric  is  understood  to  indicate  that  a  distribution  of 
loading  or  angle  of  attack  is  equal  in  absolute  magnitude  on  each  half 
of  the  wing  but  of  opposite  sign. 

^Measured  parallel  to  the  plane  of  symmetry. 
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induced  drag  coefficient 


rolling-moment  coefficient 


induced  draf 


rolling  moment  due  to  rolling  ^ ^ — r  >  radian 

Ld(pb/2V)  J 


rolling  moment  due  to  aileron  deflection 


,  per  radian 


/  PAfr  ^ 

spanwise  loading  coefficient  for  unit  rolling  moment  ^ 


scale  factor 


factors  of  loading  interpolation  function 


spanwise  loading  coefficient  or  dimensionless  circulation 

/r  \ 


^)or  r-2. 

2b  /  VbV 


spanwise  loading  coefficient  due  to  rolling ^ per  radian 

spanwise  loading  coefficient  due  to  aileron  deflection  > 

per  radian 

wing  geometry,  compressibility,  and  section  lift— ciirve— slope 


hji  integration  factors  for  spanwise  loading  due  to  ailerons 

M  Mach  number 

m  arbitrary  number  of  span  stations  defined  by  n  =  cos  S2L- 

m+1 

p  rate  of  rolling,  radians  per  second 

pb/2V  wing— tip  helix  angle,  radians 

Pvn  coefficient  depending  on  wing  geometry  and  indicating  the 

influence  of  antisymmetric  loading  at  span  station  n  on 
the  downwash  angle  at  span  station  V 
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q  free-stream  dynamic  pressure ^  pounds  per  square  foot 

S  area,  square  feet 

t  ^ratio  of  aileron  chord  to  wing  chord 

V  free— stream  velocity,  feet  per  second 

w  induced  velocity,  normal  to  the  lifting  surface,  positive  for 
downwash,  feet  per  second 

y  lateral  coordinate  measured  from  the  wing  root  perpendicular  to 
the  plane  of  symmetry,  feet 

ay  ^section  angle  of  attack  at  span  station  V,  radians 

^angle  of  antisymmetric  twist  of  the  elastic  wing  produced  by  the 
loading  due  to  rolling,  radians 

—  ^rate  of  change  of  wing-section  angle  of  attack  with  control-sur- 
face  angle  for  constant  section  lift  coefficient 

P  compressibility  parameter 

F  angle  of  sideslip,  radians 

r  dihedral  single  measured  perpendicular  to  the  plane  of  symmetry, 

radians 


spanwise  circulation,  feet  squared  per  second 

5  ^ angle  of  deflection  of  full  wing-chord  control  surface,  radians 

6  angle  of  deflection  of  full-wing-chord  control  surface,  measured 

perpendicular  to  the  hinge  line^  radians 

T]  dimensionless  lateral  coordinate 


dimensionless  aileron  span 


aileron  span 


spanwise  center  of  pressure  on  one  wing  panel 


0  trigonometric  spanwise  coordinate  9,  indicating  the  edge  of  the 
'  aileron  span,  radians 


^See  footnote  2  p.2. 
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K 


V 


A 


A 


cp 


n,V 


k 

a 

t 

T 

R 

av 


ratio  of  section  lift-curve  slope  at  a  span  station  v  to 
both  at  the  same  Mach  number 


sweep  angle  of  the  wing  quarter— chord  line^  positive  for  sweep- 
back,  degrees 


compressibility  sweep  angle  parameter 

taper  ratio  f  ^E-Chord  \ 

\  root  chord  / 


tan' 


— 1  / tan/^ 


,  degrees 


trigonometric  spanwise  coordinate  (cos“^  r\) ,  radians 


Subscripts 


integers  pertaining  to  specific  span  stations  given  by 

nit  V  jt 

T)  =  cos  —  or  T]  =  cos  ^ 
o  o 

pertaining  to  span  station  k 

center  of  pressure 

aileron 

pertaining  to  fraction-of-^ing-chord  ailerons 
wing  tip 
wing  root 
average  or  mean 


DEVELOPMENT  OE  METHOD 


The  simplified  lifting-surface  method  used  herein  replaces  a  lift¬ 
ing  sxirface  by  a  lifting  vortex  located  at  the  wing  one-quarter-chord 
line.  The  boimdary  condition  for  determining  the  vortex  strength  dis¬ 
tribution  specifies  that,  along  the  three-quarter— chord  line  of  the 
wing,  there  shall  be  no  flow  through  the  lifting  surface.  In  effect, 
this  specifies  that,  at  the  three-quarter-chord  line,  the  ratio  of  the 
velocity  normal  to  the  mean  caiilier  line  (induced  by  the  bound  and  trail¬ 
ing  vortices)  to  the  velocity  of  the  free  stream  shall  equal  the  sine 
of  the  angle  of  attack. 
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Span  loadings  are  theoretically  additive.  Since  the  symmetric 
angle-of-ettack  distribution  contributes  only  to  symmetric  loading,  it 
follows  that  the  antisymmetric  loading  is  independent  of  symmetrically 
distributed  wing  twist  or  camber;  hence,  to  find  antisymmetric  loading, 
it  is  only  necessary  to  consider  the  loading  resulting  from  the  anti¬ 
symmetric  distribution  of  the  angle  of  attack  across  the  wing  span.  In 
the  subject  case,  such  a  distribution  is  experienced  by  the  wing  as 
induced  angle  due  to  rolling^,  the  effective  twist  due  to  aileron 
deflection,  or  sideslip  of  the  wing  with  dihedral. 


For  an  antisymmetric  angle— of— attack  distribution,  the  loading  dis¬ 
tribution  will  be  equal  in  absolute  magnitude  on  each  semi  span,  but  of 
opposite  sign.  The  loading  therefore  needs  only  to  be  found  over  the 
semispan,  and,  since  the  loading  is  zero  at  the  wing  root,  only  span 
stations  outboard  need  be  considered.  The  mathematical  development  of 
the  simplified  lifting-surface  method  for  the  case  of  antisymmetric 
loading  is  given  in  Appendix  A.  As  shown  in  Appendix  A,  (in-l)/2  linear 
equations  in  terms  of  loading  distribution  are  obtained  which  satisfy 
the  wing  angle-of-attack  conditions^  at  the  three-quarter-chord  line  at 
m  stations  n,  where  m  is  an  arbitrary  odd  integer.  These  equations 
are  represented  by  the  summations 


av 


m— 1 


n=l 


Pvn 


V  =  1,  2,  3, 


mr“l 

2 


(1) 


^n  considering  the  case  of  the  angle  induced  by  rolling  as  equivalent 
to  an  antisymmetric  distribution  of  twist,  it  must  be  noted  that 
account  should  be  taken  of  the  fact  that  a  rolling  wing  leaves  a 
twisted  vortex  trail;  whereas  a  twisted  wing  does  not.  The  difference 
in  induction  effects  on  the  wing  of  the  straight  and  twisted  vortex  is 
considered  insignificant  here,  as  has  been  assumed  in  other  analyses. 
‘%ie  reader  should  note  that  the  boimdary  condition  is  given  by 
Wy  =  V  sin  from  which  (w/v)v  is  seen  to  equal  sin  ay.  The  sub¬ 
stitution  of  av  for  sin  av  has  the  effect  of  increasing  the  value 
of  loading  on  the  wing  above  that  necessary  to  satisfy  the  boundary 
condition.  However,  the  boundary  condition  was  fixed  assuming  that  the 
shed  vortices  moved  downstream  in  the  extended  chord  plane.  A  more 
realistic  picture  is  obtained  if  the  vortices  are  assumed  to  move  down¬ 
stream  in  a  horizontal  plane  from  the  wing  trailing  edge.  It  can  be 
seen  readily  that,  if  this  occurs,  the  normal  coir^onent  of  velocity 
induced  by  the  trails  at  the  three-quarter— chord  line  is  reduced  and, 
if  the  boundary  condition  is  to  continue  to  be  satisfied,  the  strength 
of  the  bound  vortex  must  increase.  It  follows  that  substitution  of  ay 
for  sin  ay  then  has  the  effect  of  accounting  for  the  bending  up  of 
the  trailing  vortices.  It  is  not  known  how  exact  the  correction  is, 
but  the  calculations  and  experimental  verification  show  it  to  be  of 
the  correct  order. 
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where 

av  antisymmetric  angle  of  attack  at  wing  station  V 

Pvn  coefficients  that  for  a  given  value  of  m  depend  on  wing 
geometry,  compressibility,  and  section  lift-curve  slope 

Gji  loading  coefficients  at  span  stations  n 


The  application  in  Appendix  A  of  the  present  report  is  with  m=7 
Since _ the  loading  at  the  midspan  station  is  known  to  be  zero,  considera- 

/m®®  stations:  n=l,2,3,  equal  to  wing  semispan 
positions  of  Ti  =  cos  (nn/8)  =  0.924;  0.707;  and  0.383.  Equation  (l) 
thus  becomes  ' 


n=l 


V  =  1,2,3 


(2) 


where  the  Integer  V 


pertains  to  span  station  q  =  cos  (vtt/S) 


To  obtain  the  loading  coefficients 
only  to  evaluate  the  coefficients  py^ 
the  antisymmetric  angle  of  attack  ay. 


^n  ~  ('^Z°/2b)j^,  it  remains 
and  the  spanwise  variation  of 


Evaluation  of  the  Coefficients 


Since  m  is  chosen,  p^  becomes  a  function  only  of  wing  geometry, 
compressibility,  and  section  lift-curve  slope.  The  effects  of  compressi¬ 
bility  and  section  lift-curve  slope  are  equivalent  to  a  change  in  wing 
plan  form  and  can  be  accounted  for  by  a  proper  adjustment  of  the  pvn 
values.  As  shown  in  Appendix  B,  py^j^  can  be  conveniently  presented  as 
a  function  of  two  parameters,  namely,  a  compressible-sweep-angle  parame¬ 
ter  defined  as  Ap=  tan— i( tan  a/S)  and  a  parameter  Hy  involving  the 
ratio  of  wing  span  to  wing  chord  and  variable  section  lift-curve  slone. 
defined  by 

(3) 

where 


ratiio  of  experimental  section  lift— curve  slope  at  span  station  v 
to  the  theoretical  value  of  2jc/|3^  both  at  the  same  Mach  number 


Compressibility  and  section  lift— curve  slope  are  discussed  in 


Appendix  B. 
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cv  vlDg  chord  at  span  station  V 

The  value  dv  is  a  scale  factor  given  hy 

dv  =  0.061  for  V  =1 

=  0.23*)-  for  V  =  2  y 

=  0.381  for  V  =3 

Eq\iation  (3)  can  be  vritten  in  alternative  form  that  gives  Hv 
terms  of  wing  geometry  parameters  that  are  more  significant;  thus 


in 


1 

^  (  Hy/ Kg^y)  (  Oy/cQ^y)  _ 

(5) 


where 

ratio  of  average  section  lift-ourve  slope  to  2jc/P  both  at 
the  same  Mach  number 

i^vAav  spanwise  distribution  of  section  lift-curve  slope  for  a 
given  Mach  number 

^vAav  spanwise  distribution  of  the  wing  chord 

(PA/k^v)  compressible  aspect  ratio  and  average  section  lift-curve- 
slope  parameter 

The  term  - - 7 - :  of  equation  (5)  gives  an  effective  aero- 

(kvAbv)  (^v/cav) 

dynamic  taper  of  a  wing-  The  distribution  of  KvAav 
Mach  number,  particularly  at  transonic  speeds  (e-g.,  due  to  spanwise 
variation  of  airfoil  section)  •  However,  since  the  distribution  con— 
•^2*ibutes  to  taper  effect,  the  loading  distribution  and  not  the  total 
loading  will  be  appreciably  affected. 

With  Hy  determined  from  equations  (3)  or  (5)  values 

of  Pviy  nine  in  all,  are  presented  in  figure  1  where  Pvn  is  given 
as  a  function  of  Hy  for  various  values  of  A^. 

For  the  case  of  straight-tapered  wings  with  arbitrary  section  lift- 
curve-slope  distribution  for  which  the  chord  distribution  is  specified 
by  taper  ratio,  evaluation  of  equation  (5)  is  given  in  figure  2  where 

(Ky/Kav)S\;  for  each  of  the  three  span  stations  is  shown  as  a  function  of 

( 3A/Kav) 
taper  ratio. 
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Evaluation  of  Antisymmetric  Angle-of-Attack  Distribution  ay 


The  antisymmetric  angle-of~attack  distributions  most  commonly 
encountered  are  those  resulting  from  rolling  wings^  aileron  deflection^ 
and  sideslip  of  wings  with  dihedral.  Evaluation  of  the  angle-of-attack 
.distributions  for  these  various  cases  is  outlined  in  the  sections  imme¬ 
diately  following. 

Rolling  wings.—  For  the  case  of  the  rigid  wing,  the  induced  veloc¬ 
ity  normal  to  the  wing  surface  is  equal  to  the  upwash  velocity  experienced 
by  the  rolling  wing.  Thus,  at  span  station  V 

av  =  p^  =  -(§)nv  (6) 

where  pb/2V  is  the  tip  helix  angle.  It  should  be  noted  that  the  rela¬ 
tion  given  by  equation  (6)  assumes  the  wing  structure  to  be  rigid  in  that 
the  distribution  of  is  completely  defined  by  the  linear  distribution 

of  helix  angle.  In  the  case  of  flexible  wings,  however,  the  expression 
for  must  be  modified  to  account  for  the  stfeamwise  angle-of-attack 

change  which  may  occur  due  to  bending  or  torsional  deflections.  In  this 
case, 

pb  \ 

—  +  (7) 

where  Zio-v  represents  the  modifying  influence  of  flexibility.  Normally, 
l!a^  is  not  considered  for  straight  wings  since  only  the  effect  of  tor¬ 
sion  (which  is  usually  small)  is  involved.  On  swept  wings,  however,  the 
effect  of  bending  can  cause  Axy  to  be  quite  large  so  that  the  ay 
distribution  may  be  affected  considerably.  Due  to  the  interaction 
existing  between  the  aerodynamic  and  structural  forces,  Aa  cannot  be 
determined  directly,  but  must  be  found  through  equations  of ^equilibriiim 
or  by  iteration.  With  the  loading  for  the  rigid  wing  provided,  however, 
the  iteration  procedure  becomes  relatively  easy  to  apply.  The  first 
approximation  of  ay  is  found  from  the  loading  of  the  rigid  wing  and 
further  refinements  of  may  be  found  utilizing  the  successive  loadings 

for  the  flexible  wing  as  determined. 

Deflected  ailerons.—  Where  the  spanwise  distribution  of  the  angle 
ay  is  to  be  considered  equivalent  to  antisymmetric  aileron  deflection, 
it  must  suffer  a  discontinuity  at  the  spanwise  end  of  the  control  surface. 
The  loading  when  such  a  discontinuity  is  present  can  be  duplicated  by  a 
proper  distribution  of  antisymmetric  twist.  In  Appendix  C,  the  anti¬ 
symmetric  twist  distribution  required  by  the  present  theory  to  give  accu¬ 
rate  span  loading  distribution  due  to  ailerons  is  found  with  the  aid  of 
zero-aspect-ratio  wing  theory  given  by  reference  6.  To  minimize  the 
confutation  involved,  it  is  convenient  to  consider  both  the  case  of  out¬ 
board  and  inboard  ailerons. 
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1.  Outboard  ailerons:  With  in=7^  three  different  aileron  spans 
can  be  conveniently  defined  for  the  outboard  ailerons.  For 
the  aileron  spans  measured  from  the  wing  tip  inboard^ 

the  antisymmetric  twist  distribution  required  per  unit  deflec¬ 
tion  of  full-wing— chord  ailerons^  coy/b,  is  given  by 


Case 

I 

II 

III 

^a 

0.169 

0.10^4 

0.805 

ai 

& 

1.003 

0.971 

0.998 

6 

• 

o 

H 

-4 

.996 

.991 

6 

.006 

.014 

.978 

(8) 


2.  Inboard  ailerons:  With  m=7^  three  different  aileron  spans 
can  be  conveniently  defined  for  the  inboard  ailerons.  For 
the  aileron  spans  qa  measured  from  the  wing  midspan  out¬ 
board,  the  antisymmetric  twist  distribution  required  per  unit 
deflection  of  full-wing-chord  ailerons,  is  given  by 


Case 

lY 

V 

VI 

■Qa 

0.556 

0.831 

1.000 

ai 

6 

o.o44 

0.013 

1.016 

02 

-.017 

.961 

.919 

3 

5 

1.087 

1.095 

1.101 

(9) 


Sideslip  of  wings  with  dihedral.—  For  calculating  the  rolling 
moment  caused  by  dihedral  angle  for  the  sideslipping  wing,  the  effect  of 
the  skewness  of  the  vortex  field  in  altering  the  effects  of  the  dihedral 
angle  will  be  assumed  to  be  small  (as  assumed  in  reference  5)-  The 
problem  then  simplifies  to  finding  the  rolling  moment  due  to  antisymmetric 
angle  of  attack  with  the  unskewed  vortex  field.  The  solution  to  this 
problem  is  the  same  as  for  the  ailerons  which  has  already  been  solved. 
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The  antisymmetric  distrihutioD  of  angle  of  attack  for  the  side¬ 
slipping  wing  with  dihedral  is  gi^  jn  by 

ov  =  p  r  (10) 


where 

<Xy  effective  angle— of~attack  distribution 

P  angle  of  sideslip  measured  positive  in  the  counterclockwise 

direction  from  the  plane  of  symmetry 

r  dihedral  angle 

The  wing  parameter  T  is  not  affected  by  compressibility.  Equation  (lO) 
is  approximate  for  small  values  of  p  and  F. 

For  unit  p  F  over  the  span  of  the  ailerons  considered, 

PF=6  (11) 

can  be  substituted  for  6  in  equations  (8)  and  (9)* 


APPLICATION  OF  METHOD 


For  the  cases  of  antisymmetric  angle— of-attack  distributions  result¬ 
ing  from  rolling,  aileron  deflection,  or  sideslip  with  dihedral,  it  is 
possible  to  present  a  set  of  simultaneous  equations  which  are  required 
for  the  solution  of  the  load  distribution  for  an  arbitrary  plan  form. 

With  the  loading  known,  integration  formulas  can  be  given  to  determine 
aerodynamic  coefficients. 

The  loading-distribution  coefficient  Gn  determined  from  the  solu¬ 
tions  of  the  simultaneous  equations,  are  functions  of  pvn  which  has 
been  shown  in  a  preceding  section  to  be  a  function  of  wing  geometry, 
compressibility,  and  section  lift-curve  slope.  The  aerodynamic  coeffi¬ 
cients  are  integrations  of  the  load  distribution  and,  therefore,  will 
also  be  a  function  of  wing  geometry,  compressibility,  and  section  lift- 
curve-slope  parameters.  Application  of  the  method  to  the  general  solu¬ 
tion  for  arbitrary  chord  distribution  is  outlined  and  solutions  are  pre¬ 
sented  for  the  case  of  straight  taper. 
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General  Solution 


AerodjTQ^Ti^ic  characteristics  due  to  rolling The  solutions  for  the 
aerodynamic  effects  due  to  the  rolling  wing  ‘will  be  found  and  loading > 
rolling  moment,  spanwise  center  of  pressure,  and  induced  drag  will  be 
obtained. 

1.  Simultaneous  loading  equations:  The  Vy-n.  values  are 
obtained  from  figure  1  and  table  I®  with  values  of  Ev 
given  by  equations  (3)  or  (5)* 


The  simultaneous  equations  (2),  for  the  rigid  and  flexible 
wing,  respectively,  become: 


where 


-0.92^  =  PiiGi  +  Pi^2  +  PIQG3 
— 0.707  =  P21O1  +  P22^2  P23^3 
"^•383  =  Psl^l  P32^2  P33^3 


On  = 


'n 


pb/2V 


y  (12) 


and 


-0>92k  +  +  P13G, 


-0.707  + 


“0 . 3  83  + 


Gi 


pb/2V 

pb/2V 

ACg 

pb/2V 


-  P21G1  +  P22*^2  1*23^3 


=  n  G  +  V  G  +p  G 
■^31  1  -*^32  2  -^33  3 


>  (13) 


where  Gr»  =  — v —  and  is  the  incremental  angle  of 

^  pb/2V 

attack  due  to  aeroelastic  effects. 


2.  Loading  distribution:  The  loading-distribution  coefficient 
is  giTOn.by  G  =  cjc/2b.  Other  forxas  of  the  loading  coeffi¬ 
cient  are  given  by  the  identities 


Values  of  Pyjj  beyond  the  scope  of  figiire  1  are  included  in  table 
For  values  of  larger  than  those  included  in  figure  1  and  table  i, 

the  Pyj^  curves  can  be  obtained  from  equation  (b8)  which  gives  the 
linear  asymptotes  of  the  p^^  fimctton. _ _ _ 
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(1^) 


G  =  i_  =  £i. 

2A  2A 

The  loading  is  knovn  to  he  zero  at  ^  =  0  and  1  and  is 
determined  at  three  intermediate  span  stations.  Values  of 
loading  at  other  span  stations  can  he  obtained  from  a  load¬ 
ing  function  derived  in  Appendix  B  or,  with  equations  (B23) 
or  (B24)  of  Appendix  B,  the  loading  can  he  found  at  span 
positions  y\  =  0.981^  O.831,  0.556,  and  0.195* 


3.  Rolling  moment:  The  damping— in— roll  derivative  for  the  solu¬ 
tions  of  equ^ions  (12)  or  (13)  is  derived  in  Appendix  B  and 
given  hy 


PCz 


+  0.707  (Gi  +  G3) 

•^av  16  \^av  / 


(15) 


h.  Spanwise  center  of  pressure;  The  equation  giving  center  of 
pressure  on  the  wing  semi  span  is  shown  in  Appendix  B  to  he 


^c.p.  “ 


pc^A 


av 


PA 

av 


(O.I63G1  +  0.24852  +  O.43OG3) 


.082  f  +  0.124/^ T— +0.215  f 

\^Z®av/4  V^l®av/fe  \*"l®av/£ 


(16) 


5.  Induced  drag:  The  induced  drag  is  derived  in  Appendix  B 
and  given  hy 


•^av 


Gi 


+  G2^  + 


Gs 


®  G2(Gi  +  Ga) 


(17) 


Aerodynamic  characteristics  due  to  aileron  deflection.—  The  solu¬ 
tions  for  the  aerodynamic  effects  due  to  ailerons  will  he  found  for  three 
different  spans  of  outboard  and  inboard  ailerons.  Cross  plots  of  these 
data  provide  curves  for  arbitrary  aileron  spans. 

1.  Simultaneous  loading  equations:  The  pyj^  values  are 
obtained  from  figure  1  and  table  I  with  values  of  Hy 
given  hy  equations  (3)  or  (5). 
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(a)  Deflected  outboard  ailerons:  The  aileron  spans  measured 
from  the  wing  tip  inboard  are  given  by  T]g^.  The  simultaneous 
solution  for  antisymmetric  spanwise  loading  due  to  deflection 
of  any  of  the  three  following  aileron  spans  can  be  obtained 
from  the  appropriate  set  of  the  following  equations: 


=  PiiGi  +  P12S2  +  PX3S3 

=  P2i5i  +  * 

^31^^  ^33^3 

(18) 

where 

(b)  Deflected  inboard  ailerons:  The  aileron  spans  measured 
from  the  wing  midspan  outboard  are  given  by  The 

simultaneous  solution  for  antisymmetric  spanwise  loading 
due  to  deflection  of  any  of  the  three  following  aileron 
spans  can  be  obtained  from  the  appropriate  set  of  the  follow¬ 
ing  equations: 


=  PiiGi  +  P12G2  +  Pis'^s 
=  P21§1  P22G2  +  P23®3 
=  Psi^l  ^32^2  +  P33^3 

(19) 


Case 

Y 

YI 

Ha 

0.556 

0.831 

1.000 

ai 

T 

0.044 

0.013 

1.016 

CX.2 

T 

-•017 

.961 

.979 

as 

6 

1.087 

1.095 

1.101 

where  Gj,  = 


Case 

I 

II 

III 

Ha 

0.169 

0.444 

0.805 

8 

1.003 

0.971 

0.998 

a2 

~ 

.017 

.996 

.991 

Os 

B 

.006 

.014 

.978 
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2.  Loading  distritiution:  The  spanwise  loading  distributions 
due  to  various  aileron  configurations  include: 

(a)  Full^ing-chord  ailerons:  The  loading  is  knovn  to  he 
zero  at  t]  =  0  and  1,  and  is  determined  at  three  inter¬ 
mediate  span  stations.  With  equation  (CI3)  and  tables  C6, 

Bl,  and  C7^  the  loading  can  be  found  at  span  stations 

T]  =  0.981^  0.831,  0.556,  and  0.195  for  each  of  the  aileron 
spans  considered.  With  these  given  points  and  the  knowledge 
that  the  slope  of  the  loading-distribution  curve  is  theo¬ 
retically  infinite  at  the  point  of  angle-of-attack  discon¬ 
tinuity  (aileron  spanwise  end),  the  loading  distribution  can 
be  faired. 

(b)  Constant  fraction  of  wing— chord  ailerons:  The  spanwise 

loading  of  constant  fraction  of  wing-chord  ailerons  is  equal 
to  the  product  of  the  loading  due  to  full-wing-chord  ailerons 
and  the  effective  change  of  angle  of  attack  with  aileron 
angle, da/d6.  The  factor  da/d6  is  a  function  of  the 
ratio  of  aileron  chord  to  wing  chord  t  =  change 

of  section  angle  of  attack  with  aileron  angle  da/d6  is 
presented  in  figure  3,  which  is  reproduced  from  figure  18  of 
reference  7. 

Although  figure  3  taken  from  reference  7  limits  the  Manh 
number  range  to  Mach  numbers  less  than  0.2,  this  limitation 
is  believed  to  be  unwarranted  since  theoi*y  indicates  that 
da/d5  is  unaffected  by  compressibility  for  the  two- 
dimensional  wing.  However,  as  indicated  in  reference  6, 
da/d6  is  strongly  affected  by  low  aspect  ratio  and  will 
change  appreciably  if  the  parameter  3A  becomes  much  less 
than  twoj  hence,  the  values  of  da/d6  from  figure  3  appear 
to  be  valid  for  pA  >2. 

(c)  Arbitrary  spanwise  distribution  of  aileron  chord:  The 
aileron  can  be  divided  into  several  spans  with  constant 
da/d6,  then  the  total  loading  is  the  sum  of  the  products  of 
the  full-wing— chord  loading  of  each  span  and  its  respective 
da/d5. 


"^In  using  da/dS  here,  it  should  be  noted  that  the  assumption  is  made 
that  the  effective  airfoil  section  is  taken  as  being  parallel  to  the 
plane  of  symmetry  and  that  the  section  approaches  a  two-dimensional 
section.  The  validity  of  this  assumption  can  be  questioned;  however, 
limited  checks  with  experiment  show  it  to  be  at  least  approximately 
correct . 
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3.  RolliTifT  moment:  Tlie  rolling  moment  can  be  found  for  the 
following  aileron  configurations: 

(a)  Will  1^-ing-chnrd-ailerons :  The  spanwise  loading  due  to 
aileron  deflection  cannot  be  integrated  with  sufficient 
accuracy  with  equation  (15)*  lii  Appendix  C,  a  similar 
integration  formula  is  developed  that  applies  to  each  given 
aileron  span.  Equation  (CIO)  and  table  C5  give 


where  for  each  of  the  cases  of  equations  (I8)  and  (19)  the 
hj^  values  are  given  by 


Case 

I 

II 

III 

IV 

V 

VI 

hi 

0.140 

0.139 

0.138 

0.146 

0.l4l 

0.l40 

h2 

• 

H 

VO 

.196 

.196 

.200 

.197 

.198 

.145 

.139 

.138 

.l40 

.139 

.i4o 

(Ta)  fraction  of  winst~chord  ailerons :  Eor  constant 

fraction  of  wing-chord  ailerons  with  aileron  angle  measured 
parallel  to  the  plane  of  symmetry,  the  aileron  effectiveness 
is  given  by 


!!!& .  ^ 


(21) 


(c)  Arbitrary  s panvise  distri'bution  of  aileron  chord:  The 
deflectioi^  of  ailerons  for  which  t  varies  spanwise  on  the 
wing  can  be  considered  as  an  eq.ui valent  wing-twist  distribu¬ 
tion.  The  effective  antisymmetric  twist  of  the  wing  is 
given  by 


(22) 
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where  da/d5  is  now  a  function  of  spanwise  position.  The 
antisyirmetric  angle— of-at tack  distribution  given  by  equation 
(22)  can  be  divided  into  spanwise  steps  of  constant  angle  of 
attack  and  the  total  rolling  moment  can  be  found  by  the 
summation  of  the  rolling  moment  due  to  each  spanwise  step. 

The  rolling  moments  of  the  spanwise  steps  are  obtained  from 
a  curve  of  rolling-^noment  coefficient  as  a 

function  of  unit  antisymmetric  angle  of  attack  from  the  wing 
root  outboard.  This* step  method  is  the  procedure  used  in 
reference  1. 

A  curve  of  /^av  as  a  function  of  unit  antisymmetric 

angle  of  attack  from  the  wing  root  outboard  can  be  obtained 
from  the  solutions  of  equation  (19)  for  the  cases  IV,  V,  and 
VI.  An  additional  point  can  be  obtained  from  the  solution 
of  case  III  of  equation  (l8),  applying  the  relations  (dis¬ 
cussed  later)  existing  between  inboard  and  outboard  ailerons. 

The  rolling  moment  due  to  the  twist  given  by 'equation  (22) 
can  be  obtained,  by  a  method  other  than  the  step  method,  from 
the  integral  given  by 


r"  da 

*^av  Jo 


(23) 


which  can  be  integrated  numerically  by  taking  the  graphical 
slopes  of  which  is  a  function  of  extent  of  unit 

ant isy333me trie  angle  of  attack  from  the  wing  root  outboard. 


Spanwise  center  of  pressure  and  induced  drag:  Spanwise  center 
of  pressure  and  induced— drag  integration  formula^  for  loading 
due  to  ailerons  are  not  given;  however,  equations  (l6)  and 
(17)  can  give  approximate  integrations  of  the  loading  to 
obtain  center  of  pressure  and  induced  drag. 

5.  Additional  considerations: 

(a)  Eelation  between  aerodynamic  characteristics  for  outboard 
and  inboard  ailerons:  The  spanwise  loading  distributions  due 
to  outboard  and  inboard  ailerons  bear  a  simple  relation  to 
each  other.  Since  loading  is  linearly  proportional  to  angle 
of  attack,  loadings  are  directly  additive.  Then,  for 
outboard  and  inboard  ailerons  with  the  spanwise  ends  of  the 
ailerons  at  the  same  span  station. 
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Glntoard  “  ^(ti  =  i)  “^outboard 

*Sl 

C7  =  Ci^  -  Cl 

^Inboard  ^(Tla=l)  “outboard  \  (24) 

^  1  —  rj 

^ inboard  ^outboard 


These  relations  do  not  apply  for  ’Hc.p.  ^Di  since 

these  characteristics  are  not  linearly  proportional  to 
loading. 

(b)  Differential  aileron  angles;  The  effect  of  a  differen¬ 

tial  between  aileron  angles  can  be  taken  into  account  by 
considering  the  of  each  wing  panel  as  one— half  the 

ant is;)mirne trie  results  of  equations  (20)^  (2l),  or  (23)*  The 
total  wing  rolling  moment  is  then  the  sum  of  the  products  of 
Ci^/2  given  by  equations  (20)^  (21)^  or  (23)  and  the  angle 
of ^deflect ion  of  each  aileron.  Although  the  total  rolling 
moment  can  be  found  by  this  procedure^  the  spanwise  loading 
distribution  can  be  found,  only  approximately  by  the  products 
of  the  antisymmetric  unit  loading  and  the  deflection 

of  each  aileron.  However^  the  loading  distribution  so  found 
will  be  quite  accurate  since  this  procedure  neglects  only 
the  small  change  due  to  the  induced  effects  of  the  differen¬ 
tially  different  opposite  wing  panels. 

(c)  Aileron  angles  measured  perpendicular  to  the  hinge  line; 
The  relationship  betvveen  aileron  angle  measured  perpendicular 
to  the  aileron  hinge  line  and  that  measured  parallel  to  the 
plane  of  symmetry  is  given  by 


tan  6 


tan  S 
cos 


(25) 


where 


sweep  angle  of  the  aileron  hinge  line 

5  angle  measured  perpendicular  to  the  hinge  line 

For  constant  fraction  of  wing— chord  ailerons  on  straight- 
tapered  wings ^  A-(-  is  given  by  ' 
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(i^)  (.6, 

where  t  is  the  fraction  of  wing— chord  aileron  measured 
from  the  wing  trailing  edge. 


Aerodynamic  characteristics  due  to  sideslip  of  wings  with  dihedral.— 
The  total  antisymmetric  loading  due  to  sideslip  can  he  considered  as  the 
sum  of  that  due  to  dihedral  angle  and  that  due  to  zero  dihedral  engle. 

For  the  unswept  wing,  the  rolling  moment  due  to  sideslip  for  zero 
dihedral  angle  is  generally  considered  negligible;  however,  for  the 
swept  wing,  this  effect  can  be  appreciable.  In  the  present  report,  only 
that  part  due  to  dihedral  angle  will  be  considered  for  the  swept  and 
nonswept  wings. 

1.  Simultaneous  loading  eq.iiations:  The  Vyn  values  are 

obtained  from  figure  1  and  table  I  with  values  of  given 

by  equations  (3)  or  (5)* 


The  simu^ianeous  equations  resulting  from  the  substitution 
of  6  =  pr  (see  equation  (11))  and  G  =  G/PP  in  equations 
(18)  and  (19)  are  applicable  in  the  determination  of  the 
effects  of  unit  outboard  or  inboard  dihedral  angle  over  the 
span  of  the  ailerons  considered. 


2.  Foiling  moment:  The  rolling  moment  due  to  various  dihedral 
angle  distributions  include: 


(a)  Constant  spanwise  dihedral  angle:  For  dihedral  angle 
constant  for  the  entire  wing  semispan,  the  loading  is  given 
by  the  solution  of  case  VI  in  equation  (19)  for  G  =  G/^T 
and  the  rolling  moment  from  equation  (20)  becomes 


- 1  =  M_  (0,l40Gi  +  O.I98G2  +  o.i4og„) 

'^av^  •'av 


(27) 


(t)  Gulled _wlng:  For  the  gulled  wing,  solutions  of  equation 
(19)  for  G  =  G/Fr  gives  the  loading,  and  the  rolling 
moment  from  equation  (20)  becomes 

pC^- 

- — I  =  (biGi  +  hsG^  +  ^363)  (28) 

av”^  av 

A  plot  of  the  results  of  cases  IV,  V,  and  VI  gives  the 
extent  of  unit  dihedral  angle  from  the  wing  root  outboard. 
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Then,  for  a  gulled  wing,  the  total  rolling  moment  equals  the 
sum  of  products  of  dihedral  angle  of  each  span  section  and 
the  rolling-moment  contribution  of  the  respective  span 
sections . 

(c)  Variable  spanwise  dihedral  angle:  If  V  varies  spanwise, 
the  rolling  moment  can  be  obtained  by  integration  as  in 
equation  (23).  The 

where - 5 - 

dT) 

part  (b)  above. 


integral  becomes 


r(Tj)  - ^ - dT] 

dT] 


(29) 


is  the  slope  of  the  curve  described  in 


Solution  for  Straight--Tapered  Wings 


Charts  of  aerodynamic  characteristics  for  straight-tapered  wings  can 
be  presented  in  terms  of  geometric,  compressibility,  and  average  section 
lift— curve— slope  parameters.  These  charts  provide  a  ready  means  of 
obtaining  data  directly. 

Aerodynsuaic  characteristics  due  to  rolling.—  The  application  of 
equation  (12)  for  a  constant  value  of  section  lift-curve  slope®  provides 
the  spanwise  loadings  at  span  stations  0.383^  0-707^  a^id  0*92^  which  are 
presented  in  figure  4  for  a  wide  range  of  plan  forms .  The  interpolation 
formula  of  equation  (B24)  will  give  values  of  loading  due  to  rolling  at 
span  stations  other  than  those  presented.  With  equation  (15),  the 
damping-in-roll  coefficients  PC2p/Kav  obtained  and  are  presented 

in  figure  5  for  a  wide  range  of  plan  forms. 

^Throughout  the  figures,  is  the  constant  spanwise  section  lift- 

curve  slope  or  the  average  of  a  small  variation.  For  large  spanwise 
variations  of  k  that  follow  the  function  given  by  equation  (Ell) 
developed  in  Appendix  B,  the  parameters  PA/k^-^-  and  X  can  be  replaced 

by  the  parameters  r - PA  , - -  and  ^  X,  respectively.  For  large 

X)/{l+X) 

spanwise  variations  of  K  that  do  not  follow  the  curve  of  equation 
(Bll),  the  simultaneous  equations  for  the  general  solution  can  be 
solved  for  arbitrary  distributions  of  k*  The  Hv  values  can  be 
obtained  from  figure  2.  _  ' 
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Aerodynamic  characteristics  due  to  aileron  deflection.—  The  appli— 
cation  of  equation  (19)^  case  IliT^of  equation  (18),  and  equation  (20) 
provide  aileron  effectiveness  in  the  coefficient  form  for 

several  aileron  spans.  In  figure  6,  P^Ig/'^av  plotted  against  extent 
of  unit  antisymmetric  angle  of  attack  from  the  wing  semi span  root 
outboard  for  a  range  of  wing  parameters. 

As  presented^  figure  6  gives  directly  the  effectiveness  of  full— 
wing— chord  inboard  ailerons  for  aileron  spans  measured  from  the  plane  of 
symmetry  outboard.  The  effectiveness  of  full— wing— chord  outboard  ailer¬ 
ons  for  aileron  spans  measured  from  the  wing  tip  inboard  is  given  by 
figure  6  directly  by  the  relations  of  equation  (24).  For  full-wing- 
chord  ailerons  located  arbitrarily  on  the  wing  semispan,  the  aileron 
effectiveness  can  be  obtained  directly  from  figure  6  as  indicated  in  the 
following  example  sketch. 


With  the  full-^ing— chord  values  given  above,  the  effectiveness  of 
constant  fraction  of  wing— chord  ailerons  or  ailerons  of  arbitrary  spanwise 
chord  distribution  can  be  found  through  use  of  equations  (2l)  or  (23)  with 
the  da/db  values  of  figure  3. 

Aerodynamic  characteristics  due  to  sideslip  of  wings  with  dihedral.— 
The  a^lication  of  equation  (19)^  case  Til  of  equation  (l6),  but  with 
b  =  pr,  and  G  =  G/pP^  and  the  use  of  equation  (28)  provides  rolling 
moments  due  to  dihedral  angle  for  the  wing  in  sideslip.  These  rolling 
moments  are  given  in  the  coefficient  form  PC^^/KavT  which  is  the  same 
function  of  t]  as  PC^^/k^v  ^ind  is  presented  with  PC^^/Kq^^  in  figure  6. 
Figure  6  with  equation  (29)  will  provide  the  rolling  moment  due  to 
sideslip  for  any  symimetric  spanwise  distribution  of  dihedral  angle. 

For  dihedral  angle  constant  spanwise,  the  rolling  moment  is  given  by 
the  value  at  t]  =  1  in  figure  6.  These  values  for  constant  spanwise 
dihedral  angle  are  presented  in  figure  7  as  a  function  of  aspect  ratio 
for  various  values  of  sweep  angle  and  taper  ratio. 
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DISCUSSION 


Effects  of  plaxL-form  parameters  on  aerodynamic  characteristics  for 
straight-tapered  wings  are  shown  hy  plots  against  the  various  parameters. 
Compressibility  is  discussed  and  formulas  given  for  a  range  of  plan  forms 
at  sonic  speeds.  Theoretical  considerations  and  experimental  comparisons 
indicate  the  order  of  reliability  of  the  present  theoretical  results. 


Straight-Tapered  Wings 


The  spanwise  loading  distribution  due  to  rolling  for  several  plan 
forms  is  presented  in  figure  8.  These  cxirves  are  the  result  of  applying 
figure  4  and  the  loading  interpolation  formula  of  Appendix  B.  The 


loading  coefficient  is  given  as 


to  inake  the  total 


loading  on  the  semispan  constant  and  thus  show  more  clearly  the  changes 
of  distribution  due  to  sweep  and  taper  ratio •  Figure  8  shows  large 
changes  in  loading  distribution  for  the  zero  tapered  wing.  The  effects 
of  sweep  are  generally  as  expected,  naunely,  that  sweepback  shifts  the 
loading  outboard. 

Effects  of  plan  form  on  the  rolling  moment  due  to  rolling  is  shown 
from  cross  plots  of  figure  5  which  are  presented  in  figures  9  and  10. 
For  higher  aspect  ratio,  figures  4,  9^  10  show  the  marked  lowering 

of  rolling  moment  due  to  sweep-  Figure  9  indicates  that  for  low  aspect 
ratio,  the  rolling  moment  becomes  essentially  independent  of  sweep  and 
taper.  The  taper  effects  on  rolling  moment  as  seen  in  figure  10  are 
small  except  for  values  of  taper  ratio  less  than  0.25. 

Typical  spanwise  loading  distributions  due  to  full-wing-chord 
aileron  deflection  are  shown  in  figure  11.  These  curves  were  faired 
with  the  aid  of  the  loading  interpolation  function  of  Appendix  C  and, 
at  the  aileron  spanwise  end,  care  was  taken  to  make  the  slope  large. 


Wing  geometry  effects  on  aileron  effectiveness  for  full-chord 
outboard  partial-span  ailerons  (with  aileron  angle  measured  parallel  to 
the  plane  of  symmetry)  are  given  in  figure  12.  The  geometry  effects  on 
PC^oAav  similar  to  those  on  the  damping— in— roll  coefficient.  Com¬ 
parison  of  figure  12(a)  with  figure  9  shows  that  approaches  the^ 

zero— aspect— ratio  value  in  the  same  manner  as  does  Figure  13  gives 

comparative  effectiveness  of  inboard  and  outboard  ailerons  for  swept 
wings.  As  sweep  increases,  the  difference  of  effectiveness  between 
inboard  and  outboard  ailerons  decreases  showing  that  inboard  ailerons 
for  highly  swept— back  wings  approach  the  effectiveness  of  outboard 
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ailerons.  Since  da/d6  becomes  large  rapidly  at  small  values  of  t 
(fig.  3),  then,  for  a  given  aileron  area,  narrow  full-span  ailerons  for 
swept-^ack  wings  may  be  more  desirable  than  larger-chord  outboard 
ailerons.  The  relative  effects  of  figures  12  and  13  apply  equally  well 
for  constant  fraction  of  chord  ailerons,  since  the  data  would  differ  only 
by  a  constant  factor  da/d5. 


Compressibility 


From  the  three-dimensional  1 inearl zed-compress lb le~f low  equation^ 
it  can  be  shown  that  the  effects  of  compressibility  will  be  properly 
taken  into  account  if  the  longitudinal  components  of  a  wing  plan  form  are 
Increased  by  the  factor  l/p.  Or,  alternatively,  if  the  linearized  com¬ 
pressible  flow  equation  be  divided  through  by  then  the  lateral  and 

vertical  components  of  a  plan  form  are  decreased  by  the  factor  p.  In 
both  cases,  the  incompressible  local  lift  is  increased  by  the  factor  1/P 
and  the  compressible  local  lift  coefficient  can  be  written  as  the 
parameter  Pc^. 

With  these  relations  known,  an  incompressible  theory  can  be  made 
into  a  compressible  theory  subject  to  the  limitations  of  the  linearized 
compressible  flow  equation.  The  geometric  parameters  of  a  wing  are 

simply  substituted  by  PA,Ap  =  tan""^  —  and  Pb.  With  local  lift 
coefficient  given  by  Pc^,  the  dimensionless  loading  becomes 
pc^c  c^c 

G  =  ~  wing- chord  distribution  remains  unaltered. 

The  sonic  speed  results  of  reference  6  can  be  used  as  a  limit  point 
in  the  present  theory  for  a  curve  of  the  variation  of  antisymmetric 
aerodynamic  characteristics  with  Mach  number.  The  following  equations 
apply  at  the  speed  of  sound  to  plan  forms  with  all  points  of  the 
trailing  edge  at  or  behind  the  upstream  line  of  maximum  wing  span: 


-kA 
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For  outboard  ailerons, 

^  sin®  0,  where  ^  ^ 

For  inboard  ailerons, 

^  (1— sin®  0),  where  =  cos  0 
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Reference  6  shows  that  aileron  effectiveness  at  the  speed  of  sound  is 
independent  of  the  chordwise  location  of  the  aileron  hinge  line^  pro¬ 
vided  the  hinge  line  remains  ahead  of  all  points  of  the  trailing  edge. 


Accuracy  of  the  Seven— Point  Solution  for  Ailerons 


The  prediction  of  aileron  effectiveness  for  given  aileron  spans 
with  wing  twist  determined  hy  zero-aspect— ratio  theory  at  only  seven  span 
points  to  satisfy  the  boundary  conditions  has  been  theoretically  shown  to 
be  stiff icient  by  comparing  results  with  the  computation  of  a  typical  3«5 
aspect  ratio,  45°  swept  wing  with  15  span  points  satisfying  the  boundary 
conditions.  The  process  of  finding  aileron  spans  for  the  15—point  method 
was  the  same  as  that  in  Appendix  C.  The  curves  showing  the  variation  of 
Ci^  with  aileron  span  for  the  7—  and  15— point  computations  were  identical. 

The  solution  for  the  angle-of-attack  distribution  that  includes  a  dis¬ 
continuity  can  be  compared  with  the  solution  for  the  continuous  aogle-of- 
attack  distribution  by  considering  an  aileron  such  that  the  angle— of— 
attack  distribution  is  equivalent  to  that  of  the  rolling  wing.  The  damp— 
ing— in— roll  coefficient  then  can  be  found  by  use  of  equation  (23)  which 
reduces  to  the  form 


and  integrating  by  parts 

This  relation  states  that  C^p  is  equal  to  the  area  between  a  curve  of 
figure  6  and  the  line  of  for  t|  =  1.  The  curves  of  figure  6  were 

foimd  by  the  simplified  lifting— surface  theory  with  antisymmetric  twist 
determined  by  zero-aspect— ratio  theory.  The  values  of  C^p  obtained  in 
this  manner  from  figure  6  were  identical  to  the  C^p  values  given  by 
siir5)lified  lifting— surface  theory  for  continuous  linear  antisymmetric- 
twist  distribution. 

As  further  theoretical  check,  the  values  of  rolling  moment  due  to 
constant  spanwise  dihedral  angle  are  obtained  from  15-point  computations 
in  reference  5  for  taper  ratio  equal  to  one,  with  which  the  present 
theory  for  the  7-point  method  is  in  exact  agreement. 


for  a  ^  ^  ^ 


pb 

2V 
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Comparison  of  Theoretical  and  Experimental  Eesnlts 


The  electro-magnetic  analogy  method  of  reference  8  provides  damping— 
in—roll  coefficients  for  an  aspect— ratio  range  of  unswept,  tapered  wings. 
The  results  of  the  present  theory  and  those  of  reference  8  are  compared 
in  figure  l4.  Except  for  the  taper  ratio  effects  on  CZp  the  comparison 
is  good.  The  rounded-wing-tip  values  of  C^p  given  by  MCA  Eep.  635 
(reference  l)  are  included  in  figure  14.  Since  rounded  wing  tips  gener¬ 
ally  give  values  of  C2p  about  6  percent  lower  than  straight  wing  tips, 
the  values  of  MCA  Eep.  635  appear  to  be  appreciably  too  high  for  lower— 
aspect— ratio  wings.  The  present  theory  and  the  theory  of  reference  8 
approach  the  value  given  by  the  zero-aspect— ratio  theory  of  reference  6 
quite  satisfactorily.  The  results  of  the  present  theory  may  be  further 
assessed  by  the  comparison  with  the  results  of  low— speed  experiment  as 
given  in  figure  15  for  the  range  of  plan  forms  presented.  For  further 
experimental  verification  of  the  accuracy  with  which  C^p  can  be  deter¬ 
mined  by  the  present  theory,  the  reader  is  referred  to  reference  5  which 
supports  the  theory  as  well  or  better  than  figure  15  of  the  present  report. 

The  loading  distributions  due  to  rolling  as  given  by  the  present 
theory  are  compared  in  figure  l6  with  low— speed  experimental  results  for 
a  range  of  swept  wings.  The  sweep  angle  seems  to  have  considerably  more 
influence  on  loading  distribution  as  given  by  experiment  than  the  theory 
indicates.  The  experimental  pressure  data,  however,  were  very  erratic 
and  no  firm  conclusion  can  be  made. 

Experimental  values  of  rolling  effectiveness  due  to  aileron  deflec¬ 
tion  are  compared  with  theoretically  predicted  values  in  a  correlation 
diagram  given  by  figure  17.  Included  are  the  results  of  a  wide  range  of 
plan  forms  which  do  not  vary  consistently  with  any  geometric  parameter  or 
aileron  configuration.  Sketches  of  the  plan  forms  and  ailerons  are  drawn 
about  the  points  of  correlation.  The  theory  makes  use  of  the  curve  of 
figure  3  giving  da/db  for  a  sealed— gap  aileron  over  a  range  of  'deflection 
of  ±10*^.  Experimental  results  for  aileron  deflections  greater  than  15^ 
measured  perpendicular  to  the  hinge  line  were  not  included.  The  correla¬ 
tion  points  of  figure  17  scatter  appreciably;  however,  the  mean  line  of 
the  points  does  approximate  the  line  of  perfect  correlation. 

Figure  17  does  not  account  for  effective  plan— form  change  due  to 
da/db.  Only  the  effectiveness  of  the  low— aspect— ratio  triangular  wing  of 
figure  17  is  exceedingly  in  error, which  is  the  result  of  neglecting  plan- 
form  change. 

The  plan— form  change  due  to  da/db  can,  in  part,  be  considered  anal¬ 
ogous  to  that  due  to  section  lift— curve— slope  change.  Thus,  the  total 
section  lift  of  a  wing,  the  chord  of  which  is  reduced  by  'da/db  and  which 
is  at  an  angle  of  attack  b,  is  equal  to  the  lift  of  the  wing-aileron 
section  for  which  the  aileron  only  is  deflected  at  the  angle  b.  This 
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change  in  plan  form^  unlike  the  section  lift— curve— slope  change  for  vhich 
the  chordwise  loading  remains  constant,  does  not  account  for  a  large 
change  in  chordwise  loading.  If  the  lifting  line  is  considered  to  he  at 
the  chordwise  center  of  pressure,  then,  for  part ial-^ing— span  ailerons, 
the  lifting  line  is  in  effect  broken  at  the  aileron  spanwise  end  and  the 
present  theory  becomes  invalid.  For  the  case  of  full-^»ring— span  ailerons, 
the  lifting  line  in  effect  remains  unbroken  and  lies  along  the  center  of 
chordwise  pressure.  For  this  case  the  wing  chord  can  be  reduced  by  da/d6 
to  account  for  plan-form  change;  however,  although  in  the  limit  of  zero 
aspect  ratio  the  results  are  the  same  as  those  of  reference  6,  this  pro¬ 
cedure  does  not  with  sufficient  accuracy  account  for  the  chordwise  loading 
shifting  aft  at  intermediate  aspect  ratios.  For  control  surfaces,  the 
effective  plan— form  change  due  to  da/dS  is  appreciable  for  the  low- 
aspect— ratio  wings  such  that  in  the  limit  of  zero  aspect  ratio  the  span- 
wise  loading  is  independent  of  the  ratio  of  aileron  chord  to  wing  chord 
(reference  6).  However,  for  moderate  aspect  ratios,  da/d6  can  be  used 
without  accounting  for  plan— form  changes  as  comparison  with  experiment 
indicates . 

Experimental  values  of  are  not  compared  with  the  present 

theory  since  reference  5  gives  ample  support  of  the  theory. 

CONCLUDING  BEMAEKS 


The  determination  of  antisymmetric  loading  for  arbitrary  wings  is 
shown  to  be  easily  obtained  by  the  solution  of  three  simultaneous  equa¬ 
tions.  The  coefficients  of  the  simiiltaneous  equations  are  presented  in 
charts  of  parameters  that  include  wing  geometry,  compressibility,  and 
section  lift-curve  slope  as  arbitrary  quantities.  Thus  the  loading  for 
an  arbitrary  antisymmetric  angle-of -attack  distribution  can  be  simply 
found  once  the  angle-of-attack  distribution  is  chosen. 

For  the  important  cases  of  antisymmetric  loading,  roll,  and  aileron 
deflection,  the  angle-of-attack  distribution  is  given  and  the  simultaneous 
equations  are  formed.  Loading  for  these  cases  can  be  found  by  simply 
obtaining  from  charts  the  coefficients  corresponding  to  the  wing  geometry, 
Mach  number,  and  lift— curve  slope,  inserting  in  the  appropriate  equations 
and  solving. 

Integration  formulas  for  the  loading  distributions  are  given  which 
enable  the  aerodynamic  coefficients  C^p  and  to  be  found.  The 

rolling  moment  due  to  sideslip  of  a  wing  with  dihedral  is  shown  to  be 
equivalent  to  that  of  aileron  deflection  and  a  procedure  for  determining 
its  value  is  given. 

For  the  special  case  of  straight-tapered  wings,  the  loading  distri¬ 
butions  and  values  of  and  are  given  in  the  chart  form  for  a 

range  of  wing  plan  forms. 
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Experimental  and  theoretical  verification  of  the  theory  is  shown  to 
he  good.  The  theory  is  applicable  for  large  aerodynamic  angles,  provided 
the  flow  remains  unseparated.  The  compressibility  considerations  are 
reliable  to  the  speed  of  sound  subject  to  the  limitations  of  the  linear¬ 
ized  compressible— flow  equation. 


Ames  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Moffett  Field,  Calif.,  Dec.  22,  194-9. 

APPENDIX  A 


EQUATIONS  FOR  TEE  DETERMINATION  OF  ANTISYMMETRIC  LOADING 


Unsymmetric  Loading 


From  NACA  TN  1476  (reference  2),  the  aerodynamic  loading  is  obtained 
by  solution  of  the  linear  simultaneous  equations 


(The  prime  indicates  the  value  for 
where 


V  =  1,  2,  .  .  .  m 

n=v  is  not  summed) 


(Al) 


(A2) 

(A3) 


byy  and  byj^  are  coefficients  independent  of  plan  form, 

M 

-1  ri'(v,o)  fno  +  L(v,M+1) 


g 


Vn  2(M+1) 


n)f. 


np 


^i=l 


(a4) 
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The  above  equations  involve  computations  over  the  entire  wing. 
However,  if  the  loading  is  assumed  to  be  symmetric  or  antisymmetric,  the 
computations  can  be  reduced  to  less  than  half  the  work.  The  case  of 
symmetric  loading  is  developed  in  reference  2  and  the  antisymmetric  case 
is  developed  in  the  following  section. 


Antisymmetric  Loading 


For  antisymmetric  loading,  the  loading  on  each  side  of  the  wing  has 
the  same  magnitude  and  distribution  but  with  opposite  sign,  or 


=  "^m+i-v^  ®n  “  “^m+i-n 
Equation  (A1)  can  then  be  written  as 


(A6) 


m-1 

2  , 

“V  =  V  )Gv  m+i^^  Gn  (^7) 

n=l 


where  the  summation  is  only  to  since  =  0  for  antisymmetric 

2 

loading . 

With  equations  (A2)  and  (A3),  equation  (A?)  becomes 


Ov 


(^bvv 


Gv  - 


nt-l 


n=l 


2  ^VQ~^V^m-fi~n 


^  ~^cv^  (^^n  ®V,m+i-n  ^ 


(a8) 
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Now,  from  equation  (A*)-) 


-1 


Svn-€v,m+i^  2(M+1) 


(^no  ^m+i-n%o^  ^ 

2 


where 


and 


L(  V^M+l)  (^njM+i  M+i^  ^ 

2 


M 

^  L(v,^l)  ( 


^n^~^m+i-ii,  |ji^ 


H=1 


m+1 


m 


^  jii  sin  >ii<Pn  cos  ^iCP|a 


^1=1 


M+1 


Prom  eqiiation  (A9)j  ^i]|i“^m+i-n,ii  defined  as 


f*  =  f  -f„^,_^  ,, 

•n|i  141  m+i-n ^  M- 


then,  using  equation  (AlO) 


m 


f* 


“  ST 


^  III  cos  (sin  ^Ii9n  -  sin 


}ii=l 


m 


2 

m+1 


^  ^1  cos  ijLiCPj^  sin  jjiiCp^  (1+  cos 


p.i=l 


and,  since  the  terms  of  the  summation  for  odd  [ix  vanish, 

m 

I 

Hi  sin  Hi^n  cos  mcp^ 


f*. 


m+1 


m 

I 


Hi=2,4,6.. .  even 


{A9) 


(AlO) 
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From  eqizatioii  (All) , 


Combining  equation  (A9)  with  (A12)  and  defining 

®*vn  ^  >in+i-n 


then 


g*vn  = 


-1 


2(M+1) 


M±1 

2 

I 

^i=0 


L(v,h)  +  l(v,m+i-^) 


f* 


n*j. 


(A12) 


(A13) 


where  for  fi  =  0  and  equal  to  half  the  values  given  by 

equation  (All)  in  order  that  the  products  can  be  fitted  into  the  sunrnia— 
tion.  With  equation  (A13)^  equation  (a8)  can  be  written  as 

JQr-1 

2  , 

cx^  =  (2Cv  +  ^  g*w)Gv  ”  V  (^^n  “  ^  S*vn)^n  (A14) 

Z_j  Cv 

n=l 


^v^m+i-v 
^V^m+i-n 

From  reference  2, 

b  sincp^ 

(cos  <Pj^  —  cos  *?v)^  L  2(m+l)  - 

which  gives  zero  values  for  ^yji  for  even  (n— v)  values.  Then,  since 
m  is  odd. 


1.2,3,... 


mr-l 


where 


“  ^VV 


Cyn  -  ^vn 


b  ,  =  0 

v,m+i-*v 


C..  =  b 


vv 


and 
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It  shoTild  be  noted  that  L(v,m.)  simplifies  somewhat  for  the  anti¬ 
symmetrically  loaded  wing  since  now  is  only  positive  in  equation 
(A5)  .  If  only  positive  values  of  rf  are  used,  then  eqiiation  (A5)  can 
be  written  as 

L*(  V,^)  =  =  L(q,n)  +  L(T|,-if) 

=  L(V,pl)  +  L(v,M+l-ti) 

In  summary,  the  foregoing  analysis  for  the  antisymmetrically  loaded 
wing  gives  _  ^ 


Pvn  <^n 


(A15) 


where 


v  =  123... 

}  }  }  *  *  *  2 


Pvtl  = 


(A16) 


=-2C,,r,  +  —  S*,m  for  n/v 


vn  ■  Cy  °  vn 


^vn  “  ^vn  ^v,m+i-n 


4  sin  cpv 

^  sin  _ 

Vn  (cos  (Pq  —  cos  <Py)^ 


l-(-l)^-v 

2(m+l) 


g*VV  =  g*\ti  ^=V 


2(M+l) 


fV 


nh  ~  m+1 


Pi  sin  picpn  cos  piCf^ 


Pi=even 
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^*no  =  f  or  =  0 
M+i  =  ^  f  or  n  = 


For  M=in,  simplifies  to 


11 

o 

O 

cp^  where 

9v 

Vtc 

m+1 

■n  =  cos 

cp  where 

- 

M+1 

=  cos 

c|^  where 

9 

nn 

n 

m+1 

For  a  discussion  of  the  relative  accuracies  obtained  for  a  choice 
of  values  of  M  and  m,  see  reference  9-  The  most  favorable  applica¬ 
tion  is  with  M=mo 


3^ 
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APPEiroiX  B 

DERIVATION  OF  RELATIONS  USED  IN  THE  METHOD 


Application  of  Appendix  A 


With  Appendix  A,  the  anti synmetri cal  loading  on  a  plan  form  for  any 
anti  symmetrical  distribution  of  ay  can  be  found.  The  principal  work 
in  the  computations  is  to  obtain  the  coefficients  of  the  simultaneous 
equations  (A15) .  These  coefficients  can  be  presented  in  charts  for  the 
complete  range  of  geometric  plan-form  parameters  into  which  are  intro¬ 
duced  the  effects  of  compressibility  and  section  lift-curve  slope.  With 
the  loading  due  to  rolling  known^  the  coefficients  and  derivatives  are 
obtained  by  integration  formulas. 


Section  lift-curve— slope  effect.—  For  a  two-dimensional  wing  with 
the  loaded  line  at  the  quarter— chord  position,  the  position  x  aft  of 
the  loaded  line  where  the  induced  downwash  equals  the  angle  of  attack  of 
the  wing  can  be  obtained  by  the  Biot  Savart  Law  as 


or 


or 


Tn  „  C7CV 

w  =  where  Tq  =  — 

2Jtx  2 

2.211. a 

V  lutx 

knx 

da  c 


then 


dcx_ 

da 


where  dc^/da  is  the  section  lift-curve  slope.  Two-dimensional  section 
compressibility  effects  that  do  not  follow  the  Prandtl-Glauert  rule  can  be 
given  consideration  by  taking  the  ratio  of  {dci /da)  at  a  given 

Mach  number  to  2jt/p.  Let  K  be  the  ratio  of  the  section  lift— curve  slope 
at  a  given  Mach  nximber  to  2it/p  or  (^c^/da)  =  2m/^y  then 

X  =  k(c/2) 

Then  the  induced  angle,  k  (c/2)  aft  of  the  loaded  line,  is  equal  to  the 
angle  of  attack  of  the  wing.  For  K  =  1,  this  is  at  the  three-quarter- 
chord  line.  For  section  lift— curve  slope  less  than  2jr,  k  is  less  than 
one  and  the  downwash  is  equal  to  the  angle  of  attack  at  some  point  between 
the  one— quarter—  and  three-quarter-chord  line. 
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To  take  into  account  the  section  lift-curve-slope  variation  in  the 
pres'^nt  theory,  the  dovnwash  must  he  found  at  a  distance  k:(c/2)  aft  of 
the  jLOj3.ded  line.  From  the  formulas  of  the  summation  in  Appendix  A, 
h/cy  should  he  taken  as  h/Kv^v^  where  Ky  is  the  ratio  of  section 
lift-curve— slope  for  a  given  Mach  number  at  span  station  v,  to  2it/|3, 

Derivation  of  parameters  for  Pyn.-—  The  Pyn  coefficients,  as 
defined  hy  equation  (A16)  in  Appendix  A,  depend  on  plan— form  geometry  in 
the  (b/cv)L*vii  functions  only,  or  pya  is  a  function  of  h/cy  and 
sweep  angle.  As  previously  shown,  h/cy  is  also  a  function  of  the  span- 
wise  variation  of  section  lift-curve  slope  and  is  effectively  equivalent 
to  h/KyCy,  where  Ky  is  the  ratio  of  section  lift— curve  slope  for  a 
given  Mach  number  at  span  station  v  to  23t/P.  The  py^^^  coefficients 
can  he  plotted  against  h/KyCy  with  sweep  angle  as  a  parameter;  however, 
h/Kycy  will  vary  from  zero  to  very  large  values  for  a  range  of  plan- 
form  geometry,  and  the  plots  become  unwieldy.  For  a  range  of  aspect 
ratio,  the  values  of  h/KyCy  are  a  maximum  for  the  zero  tapered  wings 
when  Tjy  >  0.5  (provided  plan-form  edges  are  not  concave)  and  a  maximum 
for  the  inversed— tapered  wings  for  T)y  <  0»5o  The  ratio  of  h/KyCy 
for  T]y  >  0.5  for  any  plan  form  to  those  of  the  zero  tapered  wing  or 
the  ratio  of  h/Kyoy  for  “H  y  <  0.5  for  any  plan  form  to  those  of  the 
1  nv era ed-tape red  wing  gives  a  geometric  parameter  for  any  plan  form  that 
has  maximum  values  that  depend  only  on  aspect  ratio. 

The  chord  distribution  for  straight— tapered  wings  is  given  hy 


h  A(1+A) 

cv  ‘  2  [1-  hyl(l-A)] 


(Bl) 


Then,  for  A  =  0, 


and  for  A  =  1.5 


h  _  1 

^'2(1-  hj) 

Acv  2(2+  |t)  V  I  ) 


(B2) 

(B3) 


The  ratio  of  b/K^Cy 
a  geometric  parameter  as 


to  equations  (B2)  and  (B3)  gives,  respectively. 


h/KyCy 

(b/Acv)x=o 

b/kyCy 

(1^/ACv)a=i.5 


for  0.5 


2(2+t1v) 

5 


/  b 


for 


<T)v  <  1 

0^  'H  V  <0.5 


(b4) 
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Let  Hv  “be  defined  as  two— fifths  times  the  values  of  equation  (b4) 
(the  fraction  two-fifths  is  introduced  to  give  the  approximate 

values  of  p^  to  simplify  plotting  procedures),  then  adding  effects  of 
compressibility  (see  Discussion  section) 


Hv  = 


(B5) 


where 


dv  = 


5 


for  0.^< 


<1 


^(2+t1v) 

25 


for  0 


0.5 


For  tapered  wings ^  Ey  sln^jllfies  to 


Hv 


2(1-Tiy)(l+X)  /pa's 

5[l-tiv(i->')l  \'^V^ 


for  0  •  5  —  n  V  ^  ^ 


2(2+nv)(l+X) 

25[l-riv(l-X)] 


for  0  <T]^<  0.5 


(b6) 


j 


Plots  of  pyjj^  against  Hy  in  the  range  of  Hy  =  0  to  4  will  give 
^Vh  coefficients  for  wings  of  any  chord  distribution  for  aspect  ratios 
up  to  10  or  12. 

Linear  asymptotes  of  p^yn.—  For  large  values  of  the 

functions  become  linearly  proportional  to  Ey  Since  this  linear  char¬ 
acteristic  appears  at  relatively  low  values  of  Hy,  the  simply  linear 
relation  between  py^^  and  Ey  is  quite  usable. 
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The  function  of  Appendix  A  Is  mialtlplled  hy  b/cv  and  the 

product  is  linearized. 


(  L*  =  2  ^  ^ 

\  Oy  /  cos  A  \p^  / 


-  ...ai . - 

2  -2 
^  -Tl 


+ 


sin  A- 


2ti 


tan  A 


for  T|  <  “q 


1  1  \  .  1  ^  A 

-  +  —  1  sin  A - tan - 

h-^l  n  ^  2t|  2 


...  for  q=q 

2  ^  ^  2  sin  A 

cos  A  \cv/  ^  q 


for  q  =  0 


> 


with  the  values  of  equation  (B?)  substituted  Into  equation  (A16) 
the  values  of  py^  for  arbitrary  sweep  angle  are  obtained.  Thus, 
for  m=7,  the  following  equation  (b8)  gives  values  for  p^^^^  as 
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+  1.026  tan  a)  Hi  +  7.968  -  1.494  sin  A  +  0.0l4  tan  ^  + 
COS  A  / 

0.082  (  tan^._A^  _  (l- V  1  ^  0».0.1T7_t^ 

V  tan  K  ^  ^  tan  A  ^ 


0.03^ 


1  +  0,1717  tan^A 
tan  A 


-  2.901  tan  a)  Hi  -  3.138  +  1.080  sin  A  -  0.03^  tan  ^  - 

COS  A  / 

0.034  +  0-0016  tan^  A  y_n  A-  A  +  0.1717  Wa' 

\  tan  A  '  ^  A  / 

+  1.026  tan  aJ  Hi  +  0.129  -  0.869  sin  A  +  O.082  tan  ^  + 

0  014  A--/  0-0016  tan^A\  ^  A-  Vl  -f  0.0177  tanf a)  H- 

V  tan  A  /  ^  tan  A 


0.034 


0.221 

cos  A 


0.088 


0.037 


1— V  1  +  0.1717  tan^A 
tan  A 


+  0.53^  tan  a)H2  -  2.088  —  0.383  sin  A—  0.018  tan 


1-71  +  0.0177  tan^A 


tan  A 


-  O.Okk 


VI  +  0.029^  tan^A\  _ 
tan  A  / 


V  1  +  0.0886  tan^A 
tan  A 


H2  +  4.596  -  0.146  sin  Ar-0.o44  tan  5-  + 
cos  A  ^ 


0.125 


0.125 


*1—  a/~T  +  0V0177  tan^A\  ^  Q  ^ 1““  +  0.0294  tan^A  j  _ 


tan  A 


tan  A 


i—  /y/  1  +  0.0886  tan^A'^ 


tan  A 
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0.221 
COS  A 


0.534  tan  a)H2  -  1.912  +  0.221  sin  A  +  0.10?  tan  — - 
/  2 


0.044 


1-V  1  +  0.0177  tan^A 
,  tan  A 


+  0.018 


1—  V  1  +  0.029^  tan^  A 
tan  A 


o.o44 


1-y  1  +  0.0886  tan^  A 
tan  A 


=  ( =2i22§  _  0.164  tan  a)  H,  +  0.l49  +  0.324  sin  A  +  0.034  tan  ^ 
\cos  h  /  ^  2 


0.082 


,1  +  O.I7IT  tan^A  >,  fk 

tan  A  )  \ 


-  V  1  +  0.0886  tang  A 
tan  A 


0.197 


1  +  0.1993  tang A 

tan  A 


_  /O.136 
Vcos  A 


+  0.464  tan  A )  Ho  -  1.570  -  O.389  sin  A  -O.082  tan  -  + 


1  +  0.1717  tan^A  \  Q  A-  A  +  0.1993  tan^A 

V  tan  A  J  \  tan  A 


'0.628 

cos  A 


0.082 


0.03^ 


—  0.164  tanA  1  Ho  +  3.4l7  +  O.O83  sin  A +  0.197  tan  — 


1-/7  1  +  0.1717  tan^A  V  0  163  1-  7  l  +  0.0886  tan^~A  ^ 

tan  A  y  *  V  tan  A  / 


1-/7  1  +  0.1993  tan^  A 
tan  A 


Linear  spanwise  distribution  of  ( Kc)y /( Kc)Q^.y..—  With  the  condition 
that  ih.e  ^oduct  of  section  lift-curve  slope  anowing  chord  varies 
linearly  spanwise^  then 


[1+(k^/kj^ )X  ] 


l-Tlv[l-{Kj/Kg  )X,  ] 


ho 
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and  equation  (3)  becomes 


l+{  Kfji/ 


2  -T  ^ 


{B9) 


vhere  A^  is  the  aspect  ratio  based  on  the  ving  chord  equal  to  Kc. 

In  equation  (B9),  Hv  is  reduced  to  terms  of  two  parameters.  Expressions 
of  A^  in  terms  of  aspect  ratio  for  straight— tapered  wings  and  the  dis¬ 
tribution  of  section  lift-curve  slope  can  be  found. 


Eor  straight— tapered  wings 


2b 

( 1+^, ) 


and  since  kc  is  linear 


2b _ 


then 


(Kj^+k^)/1+X 


and  equation  (B9)  becomes 


Hy  -  dy 


PA 


l+( Kj/ Kg)X 


(BIO) 


The  distribution  of  k  for  straight— tapered  wings  is  given  by 


'V 


(kc), 


[1— (  ] 

(l— x<) 


(Bll) 


Equation  (BIO)  is  in  terms  of  two  parameters  given  by  !  - - — - - 

L  (Kj^-WrjtX)/l+X  _ 

and  (k)p/kj^)X.  Solutions  for  spanwise  loading  in  terms  of  these  two 
parameters  and  Ap  are  valid  for  the  distributions  of  section  lift- 
curve  slope  given  by  equation  (Bll).  Equation  (Bll)  indicates  that  at 
X  =  ly  K\}  is  a  linear  function  and  at  X  =  0,  Ky  is  a  constant. 


NACA  TN  2140 


kl 


For  values  of  X  between  0  and  1,  Ky  is  a  curve  in  the  region 
between  the  linear  function  and  a  constant. 


Equation  (BIO)  is  given  by  figure  2  for  m=T,  but  with  the  ordi- 

Hv 


nate  given  by  the  parameter 


Pa/[(Kj^+k^X)/(1+X)] 


and  the  abscissa  by 


For  the  case  of  linear  distribution  of  (kc)v  and  straight- 
tapered  wings  for  which  the  chord  and  section  lift-curve  slope  can  be 
specified  in  three  parameters^  the  loading  and  associated  aerodynamic 
characteristics  can  be  presented  for  a  range  of  the  parameters  Ap, 

(3a/[  (kj^+KqiX)/(1+X)  ]  ,  and  (Kip/Kj^)Xo 


Integration  of  Antisymmetric  Loading 


Rolling-moment  coefficient  and  derivatives.— 
ficient  is  given  by 


F  G(n)n<iTi 


Rolling-moment  coef— 


(B12) 


where 

T\  =  COS  9 

which,  by  an  integration  formula, 


m 


m+1 


=  ITT  7  ^(^n)  sin 


n=l 


(B13) 


becomes 


pc,  = 


«PA 


2(m+l) 


n=l 


G  cos  CP  sin  cp 
n  ^n  ^n 


m 


/  G  sin  2cp 
4(m+l)  ^4i 
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Since  the  loading  is  antisymmetric,  0,  and 


mr-l 

2 


PCi  = 


«PA 


2(m+l) 


Gn  2cp^ 


n=l 


(Bl4) 


For  spanwise  loading  due  to  rolling,  the  loading  is  found  as  a 
function  of  pb/2V;  then  equation  (B14)  divided  by  pb/2V  gives 


PCZp 


itpA 

2(m+l) 


m—l 

2 


^  Gn  sin 


n=l 


(B15) 


where  _ 

G  =  G/(pb/2V) 


The  rolling  moment  due  to  ailerons  will  be  found  in  Appendix  C. 

Induced  drago—  The  induced  drag  coefficient  is^  with  equation  (B13)^ 
given  by 


pCp.  =  PA  o-iGdri  =  ^  ^  Gv'Hv  '■'v 

v=l 


where  aj^  is  one— half  the  induced  angle  of  the  wing  wake  given  by 
equation  (Al4)  for  Cy  =  then  for  antisymmetric  loading 


27tPA 

m+1 


m-1 

2 


V=1 


^VV^V 


mr-1 


n=l 


sin  cpy 


(B16) 


where  the  prime  indicates  the  value 

Spanwise  center  of  pres sure >— 
half  panel  is  given  by 


of  n  =  V  is  not  summed. 

The  center  of  pressure  on  the  wing 

GqdT] 
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The  nxmerator  is  equal  to  |3C^/|3A.  If  the  Fourier  series  for  loading  is 
assumed. 


G(cp) 


m 

=  I 

|ii=even 


sin  |ii  cp,  the  denominator  becomes 


m 

I 

|ii=even 


Vi 


sin  cp  sin  M-icpdcp 


I 


|ai=even 


then 


|ai=even 


(BIT) 


where  a^j_  are  the  Foin’ier  coefficients • 

Load ing-Hiue— to-rolling  function  and  interpolation  table.—  The 
Fourier  series  that  approximates  the  antisymmetric  loading  with  only  a 
few  terms  is  given  by 


G(cp) 


)-Li=even 


sin  piCp 


(B18) 


The  loading  is  determined  at  span  positions  of  ti  =  cos  cp^  where 

nit  ^ 

^n  =  The  a^^  are  given  by 


a|ii 


sin  M-iCpdcp 


(B19: 


With  the  quadrature  formula  of  equation  (BI3),  equation  (BI9)  becomes, 
for  antisymmetric  loading, 


nt-l 

k  V" 

n=l 


(B2( 
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For  m  =  7,  the  coefficients  are  equal  to  (for  even  (ii) 

a,  =  i  (dl.  Gi  +  Ga  +  —  Os') 


2  V2 


a^  =  —  (Gi  —  G3) 


(B21) 


^6  =  1  (4^  Gi  -  Ga  +  -^  Gg^ 

Equation  (BI8)  with  (B2l)  can  be  arranged  to  give 

G  (cp)  =  sin  2  cp  +  sin  4  cp  +  dk.  sin  6  9 j  Gi  + 

2\  2  2  / 

i-  ( sin  2  cp  -  sin  6  cp  )Ga  + 

sin  2  cp  -  sin  4  cp  +  dk  sin  6  cp^Gg 


(B22) 


With  equation  (B22)  the  loading  due  to  rolling  can  he  determined  at  any 
span  position.  Letting  cp  =  cpj^  =  ^  and  tabulating  the  factors  of 

©lik^  ^  interpolation  table  may  be  obtained  to  determine  loading  at 
span  station  k, 

TABLE  Bl, 

[iri=T] 


where 


nk 

0.981 

0.831 

0.556 

0.195 

n^\^ 

1/2 

3/2 

5/2 

7/2 

1 

0.8155 

0.5449 

-0.1622 

0.1084 

2 

-.2706 

.6533 

.6533 

-.2706 

3 

.1084 

-.1622 

.5449 

.  .8155 

G]j  =  ^  e-nkCr-n 

n=l 


(B23) 
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Equation  (B23)  may  be  used  for  interpolation  of  any  form  of  loading  coef- 
ficient,  thus 


( 


\ 

k 


’nk 


c^c 


n 


(B24) 


APBEKDIX  C 

DETERMINATIOIT  OF  ANTISYMMETRIC  WING  TWIST  FOE  FINDING 
SPANWISE  LOADING  DUE  TO  AIIERON  DEFLECTION 


Wing  T^'^ist  for  a  Given  Aileron  Span 


The  determination  of  loading  for  an  angle— of— attack  distribution 
that  contains  a  discontinuity  by  a  method  which  satisfies  the  boundary 
conditions  at  a  finite  number  of  points  can  be  made  by  increasing  the 
number  of  points  until  the  solutions  become  sufficiently  accurate «  For 
the  method  as  given  in  Appendix  A,  the  number  of  points  that  satisfy  the 
boundary  conditions  is  given  by  mo  For  the  large  value  of  m  required 
for  accurate  results^  the  computations  become  exceedingly  laborious;  how¬ 
ever^  a  procedure  using  a  moderate  value  of  in  can  be  determined  by  use 
of  a  low-aspect— ratio  theory  with  which  a  wing  twist  can  be  found  that 
duplicates  the  results  of  the  discontinuous  angle-of-attack  distribution. 

A  theoretical  but  relatively  simple  method  of  finding  spanwise  load¬ 
ing  due  to  inboard  and  outboard  ailerons  for  wings  of  low  aspect  ratio  is 
given  by  reference  6.  In  the  present  theory,  as  aspect  ratio  approaches 
zero,  g* Va  values  of  Appendix  A  become  zero  and  the  Pvn  coefficients 
given  by  equation  (A16)  become  constant  or  independent  of  plan— form  shape 
and  equal  to 


Pvn  =  -2Cvn 
Pvv  “ 


(Cl) 


These  coefficients  are  given  by  the  relations  under  equation  (A15)  and 
Pvn  tabulated. 
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TABIjE  Clo“  P  Vn 
[For  m=7  and  A=0] 


n 

1 

, 

2 

3 

1 

10.4524 

-2.0000 

0 

2 

-3.6954 

5.6568 

-1.5308 

3 

0 

-2.0000 

4.3296 

With  equation  (A15)j  antisymmetric  loading  can  he  found 'for  zero— aspect- 
ratio  wings.  As  a  comment  on  the  accuracy  of  the  present  theory  for 
m=7,  the  solution  of  equation  (A15)>  with  the  A=0  py^  values  for  load¬ 
ing  due  to  rolling  gave  the  same  values  at  the  three  semi span  stations  as 

(pb/2V)  sin  2cp 

does  reference  6,  namely^  G(cp)  - - - . 

The  zero-aspect— ratio  theory  of  reference  6  shows  that  all  span 
loading  characteristics  are  independent  of  plan— form  shape  for  zero 
aspect  ratio.  This  independence  makes  that  theory  ideal  for  obtaining 
the  boundary  conditions  of  the  present  theory  for  zero  aspect  ratio, which 
should  apply  with  the  present  theory  for  higher,  aspect  ratios  for  which 
plan-form  shape  has  an  effect  on  spanwise  loading.  The  boundary  condi¬ 
tions  of  the  present  theory  are  given  by  the  antisymmetric  values  of 
in  equation  (A15).  The  problem  is  to  find  what  antisymmetrical  distribu¬ 
tion  of  ay  is  required  for  the  present  theory  to  duplicate  the  exact 
loading  distribution  given  by  reference  6  for  a  given  aileron  span. 

The  aileron  spans  are  arbitrarily  chosen  for  the  present  theory  as 
the  mean  value  of  the  spanwise  trigonometric  coordinate  of  the  downwash 
point  at  a  section  angle  of  attack  equal  to  zero.  For  m=7^  three 
aileron  spans  can  be  defined  for  both  outboard  and  inboard  ailerons. 

Let  Tia  aileron  span,  eind  9  the  spanwise  point  of  the  end  of 

the  aileron,  then 


n  =  1  —  cos  6  for  outboard  ailerons 

*8.  I 

Ti  =  COS  0  for  inboard  ailerons 

'a 
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For  the  present  theory,  the  aileron  spans  defined  are  tabulated  as  fol- 
lovs : 


TABLE  C2 


Outboard 

Inboard 

Case 

I 

II 

III 

IV 

V 

VI 

0 

3tt 

'll 

3)t 

16 

16 

16 

16 

16 

0 

^a 

0.1685 

0.4444 

0.8o49 

0.5556 

0.8315 

1.0000 

For  the  aileron  spans  listed  in  table  C2,  tlie  exact  span  loading 
distribution  can  be  found  from  reference  6.  With  the  values  listed 

in  table  Cl  and  the  exact  values  of  Gi,  and  Ga  from  reference  6, 
equation  (A15)  gives  the  twist  required  for  the  present  theory  to  give 
the  loading  distribution  for  each  case  listed  in  table  C2  or 


r-io-'*5='‘(r)-3-«5'*(r) 

^  ^  5.6568  2 

-1.5308 4.3296(^j2^ 

o  \^o  /  \  o  / 


}  (C2) 


The  spanwise  loading  distribution  from  reference  6  for  outboard  ailerons 
is  given  by 


G(cp)  ~ 

S 

outboard 


(cos  cp  —  cos  9) 


In 


In 


sin 


sin 


cos 


0+9 

2 

0-9 

2 

0+9 

2 


0-9 
cos  — 


2 


( cos  cp  +  cos  0) 


(C3) 
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For  the  full-wing-span  aileron,  6 


jr 

2  ''''  = 


p 

—  cos  CD  In 

Tt  ^ 


1  +  sin  cp 
cos  cp 


(C4) 


For  inboard  ailerons,  with  the  same  value  of  0 


'G^  ' 

_  rG(cp)i 

rG(<p)i 

6 

B 

inboard 

na=i 

6 

outboard 


(C5) 


With  equations  (C3),  (C4),  and  (C5),  the  spanwise  loading  Gi,  G2, 
and  G3  at  span  stations  cp  =  Jt/8,  )t/4,  and  3^/Q}  or  n  =  0.9239> 
0.7071,  and  0.3827  can  he  tahcilated  for  each  of  the  cases  given  in  table 
02. 


TABLE  03 


5 

Case  I 

II 

III 

IV 

V 

VI 

B 

0.1136 

0.1919 

0.2316 

0.0454 

0.1237 

0.2373 

6 

.0500 

.2800 

.3851 

.1164 

.3464 

.3964 

^3 

6 

.0190 

.1022 

.3620 

.2922 

.375^^ 

.39^ 

The  twist  distribution  required  for  each  case  is  obtained  with 
equation  (02)  and  table  03,  tabulating 

TABLE  04 


With  the  twist  distribution  given  by  table  04,  equation  (A15)  can  be  used 
to  solve  for  spanwise  loading  due  to  ailerons  for  any  of  the  six  cases. 
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Rolling  Moment  Due  to  Aileron  Deflection 


The  rolling  moment  is  given  by 

It 

G(cp)  sin  2  cpdcp  (C6) 


For  span  loading  due  to  ailerons,  the  loading  distribution  is  distorted 
sufficiently  such  that  the  quadrature  formula  given  by  equation  (B13)  is 
not  sufficiently  accurate  for  m  =  7  to  integrate  equation  (C6).  With 
equation  (B18) 

Cz  =  02  (C7) 


Expanding  equation  (B18)  for  cp  =  jt/8,  ix/k,  and  3«/8,  or  obtaining 
^ly  and.  G3  in  series  of  a*s,  the  sum  of  the  G*s  gives 

a2  =  i  (O.707IG1  +  G2  +  0.7071Gq)  +  a-L4  -  a^g  +  ago  -  a34  (c8) 

The  higher  harmonic  coefficients  can  he  put  as  factors  of  the  G^.  The 
rolling'-moment  coefficient  becomes 


A  ^ 

»  1  V  > 

16 

It 

1  + 

16 

0.7071n 

(®’14  ^18  ®^3o  Q’34) 


16 


1  + 


18 


30 


I.2O7IG1  J 

G2  + 


Gi 


I.207IG2 

(ai4  ~  a  10  4-  ago  —  334) 
1.2071G^ 


G 


When  hjx  ±&  defined  as  the  coefficients  of  G^. 

=  A  (hiGi  +  h2G2  +  hgGs) 


(C9) 


(CIO) 


The  ratio  of  (a^4  —  aig  +  ago  —  a34)  to  G^  can  be  evaluated  by 
the  zero-aspect-ratio  theory.  It  is  expected  this  ratio  vill  not  vary 
appreciably  "vrith  aspect  ratio.  From  reference  6^  the  loading  in  series 
expansion  gives  for  equation  (B18) 
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( 


^i\ 

^  '  outboard 

) 

^  S  4.=  1 


lil. 

-(-1)  ^ 


(cos  6  sin  Hi  0  -  M-I  sin  0  cos  Hi  0) 


V 

(—) 

^  ^  '^inboard 


2-1 


>(C11) 


^  ^  ^Tjg^  =r  1  ^  ^outboard 


These  high  harmonic  coefficients  are  small^  but  are  not  negligible  for 
loading  due  to  ailerons.  The  hn  are  tabulated  for  each  of  the  cases. 


TABLE  C5 


hn 

Case  I 

II 

III 

IV 

V 

VI 

111 

0.1398 

0.1388 

0.1379 

0.1462 

0.i407 

0.1402 

h2 

.1994 

.1963 

.1955 

.2004 

.1973 

.1975 

hg 

.1446 

.1388 

.1382 

.i4oo 

.1394 

.1397 

Spanvise  Loading  Distribution 


The  spanvise  loading  distributions  due  to  the  twist  distributions 
of  table  C4  are  found  at  three  span  stations,  and,  since  these  loadings 
are  not  cou^letely  defined  by  a  few  terms  of  the  assiuned  loading  series, 
the  values  of  loading  at  other  span  stations  cannot  be  found  accurately 
by  direct  use  of  equation  (B23)  and  table  Bl.  For  zero— aspect-ratio 
wings,  the  spanwise  loading  distribution  due  to  aileron  deflection  is 
given  at  all  span  stations  by  equation  (C3)«  The  loading  distribution 
for  other  than  zero-aspect— ratio  wings  will  fluctuate  about  the  value 
given  by  equation  (C3)  in  a  manner  similar  to  the  manner  that  loading 
due  to  rolling  varies  about  the  function  sin  2  cp  of  zero-aspect— ratio 
theory.  Since  the  interpolation  table  of  equation  (B23)  applies  only  to 
loadings  that  vary  about  the  function  sin  2  cp  the  loading  due  to 
aileron  deflection  can  be  divided  by  the  ratio  of  equation  (C3)  to  sin  2  cp 
and  the  resulting  loading  will  be  approximately  given  by  sin  2Cp. 
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The  zero-aspect— ratio  values  of  equation  (C3)  can  be  tabiilated  as 

ratios  of  .  Define 

sin  2cp 


En  = 


G(cto)/S 
sin  2Cpj^ 


(C12) 


The  zero— aspect-ratio  values  of  En  can  be  tabulated  for  each  aileron- 
span  case  considered. 


TABLE  C6.-  En 


Outboard 

Inboard 

Case 

I 

II 

III 

IV 

V 

VI 

'\^na 

n 

0.1685 

0.4414 

0.8049 

0.5556 

0.8315 

1.0000 

1 

0.1607 

0.2714 

0.3275 

0.0642 

0.1749 

0.3358 

2 

.0500 

.2800 

.3851 

.1164 

.3464 

.3964 

3 

.0269 

.1445 

.5119 

.4132 

.5309 

.5578 

The  interpolation  series  of  equation  (B23)  becomes 


(C13) 


where  G  =  r  and  e^^  are  given  by  table  Bl.  With  Rj,  tabulated, 

^  kTC 

values  of  loading  at  span  stations  Tjjj  =  cos  g—  are  obtained. 


TABLE  C7.-  Ei 


^NOase 

I 

II 

III 

IV 

V 

VI 

0.981 

1/2 

0.1738 

0.2663 

0.3162 

0.0559 

0.1484 

0.3222 

.831 

3/2 

.1056 

.2787 

.3499 

.0802 

.2533 

.3589 

.556 

5/2 

.0341 

.2250 

.4365 

.2424 

.4225 

.4566 

.195 

7/2 

.0233 

.1212 

.5304 

.6363 

.7343 

.7575 
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TABLE  I.-  ANTISYMMETRIC  INFLUENCE  COEFFICIENTS, 
BEYOND  THE  SCOIE  OF  FIGURE  1 


NACA  W  211)-0 


(a)  I,  n  •  ! ,  df=0.06/. 

^ure  /.-  Influence  coefficients,  ,  for  antisymmetric  spanwise  loading  plotted  as  a 
function  of  the  wing  geometric  parameter,  for  values  of  the  compressible  sweep 
parameter,  degrees. 
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Figure  /-  Continued. 
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Figure  /-  Continued. 
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0.381. 


Figure  3-  Variation  of  iift-effectiveness  parameter  with  aileron  chord  ratio, 
t  = .  Average  trailing- edge  angle  about  iO’  -,  M^O.2.  Curves  from 
reference  7 . 
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(a)  ff  =0.3827 . 


CfC 

Figure  4 -  Variation  of  loading  due-fo -rolling  coefficient 

pressible  sweep  parameter  A^,degrees,  for  straight- tapered  wings. 
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-50  -40  -30  -20  -10  O  10  20  30  40  50  60  70 


(c)  ^=0.9239. 
Figure  4  -  Concluded. 
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(a)  \=0. 


Figure  5. -Variation  of  damping-in-roii  parameter  ■  with  compressible 
sweep  parameter  A^^degrees,  for  straight -tapered  wings. 
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(b)  A  =0.25. 


Figure  5.  -  Continued. 


NACA  TN  214-0 


-50  -40  -30  -20  -to  0  10  20  30  40  50  60  70  80  90 


(c)  X  =  0.50. 

Figure  5.-  Continued. 
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(d)  A  =  1.0. 


Figure  5.-  Continued. 
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Figure  5 r  Concluded. 
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V 


(a)  <4  -  O. 


fiC2 

Figure  6 -  Aileron  rolling-moment  parameter  —;^,per  radian,  and  rolling 

0C2  " 

moment  due  to  sideslip  with  dihedral  — j:  radian  sguared,  for 

extent  of  unit  antisymmetric  angle  of  attack  from  the  wing  root 
outboard. 
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Figure  6  -  Continued 
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Figure  6.  -  Concluded 
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NACA  33f  2140 


Figure  8-  Spanwise  loading  due  to  rolling  of  wings  with  various  taper 
ratio  and  sweep  angle  parameters  of  aspect  ratio  parameter  *  4 . 
The  curve  for  A»0  serves  as  a  basis  for  comparison  and  the  fach>r 


%p 

\p(A  ’O) 


gives  the  curves  constant  area  equal  to  3.454 . 
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Figure  8.  -  Concluded. 
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NACA  IN  2ll|-0 


Figure  9-  Variation  of  the  damping-in- roll  parameter  — with  aspect 
ratio  parameter  for  various  sweep  angles  and  taper  ratios. 


Figure lO.-Vanation  of  the  damping-in- roll  parameter  - ^  with  taper 

^av 

ratio  for  wings  with  various  aspect  ratio  and  sweep  angle  parameters 
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NACA  TN  2ll«D 


Figure  U-  Antisymmetric  spanwise  loading  distribution  due  to 
deflection  of  outboard  ailerons,  for  Ca=c,  on  wings  hav¬ 
ing  a  taper  ratio  of  0.5. 


Figure  12.-  Aileron  rolling-moment-coefficient  parameter  — 

radian,  for  outboard  ailerons  as  a  function  of  aspect  ratio 
parameter,  taper  ratio,  and  compressible  sweep  parameter. 
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Igure  13  -  Aileron  rolling— moment  parameter  — 17- ,  per  radian,  for  =  A  -0.5,  as 

Kay  Kav 

a  function  of  outboard  and  inboard  aileron  spans. 


Figure  14. -  Comparison  of  damping- in- roll  coefficient  of  the 
present  theory  with  those  of  the  theories  of  references  /,  6,  ond  8. 
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Experimental 


Figure  15.- Correlation  of  theoretical  and  low- speed  experimental  damping- 
in- roll  coefficient  for  various  plan  forms. 


Figure  iS.-Comparisons  of  theoretical  and  tow-speed  experimental  spanwise 
loading  coefficients  ■  due  to  rolling  of  various  swept  wings. 
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Figure  17.- Correlation  of  theoretical  and  iow-speed  experimental  aiieron- 
effectiveness  Ci^  ^per  radian,  due  to  two  antisymmetrically  deflected 
ailerons,  for  various  plan  forms. 
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